Discrete breathers for understanding reconstructive 
mineral processes at low temperatures 

JFR Archilla* J Cuevas 
Grupo de Fisica No Lineal. Universidad de Sevilla. 
Departamento de Fisica Aplicada I. ETSI Informatica 
Avda. Reina Mercedes, s/n. 41012-Sevilla, Spain 

MD Alba, M Naranjo and JM Trillo 
Departamento de Quimica Inorganica. Universidad de Sevilla. 
Instituto de Ciencia de Materiales de Sevilla. 
Consejo Superior de Investigaciones Cientificas. 
RO. Box 874, 41080-Sevilla (Spain) 

September 20, 2006 



Abstract 

Reconstructive transformations in layered silicates need a high tem- 
perature in order to be observed. However, very recently, some systems 
have been found where transformation can be studied at temperatures 
600° C below the lowest experimental results previously reported, in- 
cluding sol-gel methods. We explore the possible relation with the exis- 
tence of intrinsic localized modes, known as discrete breathers. We con- 
struct a model for nonlinear vibrations within the cation layer, obtain 
their parameters and calculate them numerically, obtaining their ener- 
gies. Their statistics shows that although there are far less breathers 
than phonons, there are much more above the activation energy, being 
therefore a good candidate to explain the reconstructive transforma- 
tions at low temperature. 
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1 Introduction 



During the last decade, some of the present authors have achieved the syn- 
thesis of crystalline high-temperature polymorphs of rare earth (RE) dis- 
ilicates (RE2Si207) at non-expected temperatures, significantly lower than 
those previously reported, through a reconstructive process (LTRT) from 
clay minerals as the silicon source. 1-4 

This finding is of general importance in the development of advanced 
structural ceramics 5 or the storage of radioactive wastes 3 and should allow 
completion of the available Si O2-RE2O3 phase diagram. 6 Although no 
precise explanation has been found up to now, some of the present authors 
had previously suggested a chimie douce mechanism based on the diffusion 
of RE ions into the interlayer space of the expandable clay minerals. 3 ' 4 

MacKay and Aubry 7 have suggested that a possible effect of the ex- 
istence of localized nonlinear vibrations, named discrete breathers (DBs) 
could be an apparent violation of Arrhenius' law, i.e., the phenomenon of 
chemical reactions taking place at much lower temperatures than expected. 
Although this hypothesis is adventurous, it is worth exploring. Moreover, 
experimental evidence of DBs has already been found in several systems such 
as antiferromagnets, 8 waveguide arrays, 9 molecular crystals 10 or Josephson- 
junctions. 11 Moving DBs have been proposed as an explanation of dark 
tracks in muscovite 12 and there existence in muscovite has just been proven 
through a sputtering experiment. 13 

With this aim, we have made calculations and shown that the contribu- 
tion of DBs can provide an interpretation for LTRT in clay minerals. And in 
order to give experimental support to the hypothesis of DBs, it will be shown 
that the LTRT phenomenon is not exclusive of expandable clay minerals, as 
expected by the previously suggested "chimie douce" mechanism, 3 ' 4 but also 
extensible to non-expandable layered silicates, such as mica muscovite. The 
layout of this work consists of: Section 2: Some structural considerations 
on the reconstructive nature of transformation of layered silicates into disil- 
icates; Section 3: A report on a new experiment on LTRT performed by the 
authors on mica muscovite; Section 4: Argumentation about the difficulty 
to explain by the conventional chemical kinetics model the latter experi- 
ment; Section 5: Description of an alternative model based on DBs with 
numerical calculations; Section 6: After a summary of breather statistics 
theory, the description of our numerical simulations and the consequences 
on the reaction rate, we conclude with the possibility of explaining, for the 
first time, the LTRT phenomenon in the synthesis of high-temperature poly- 
morphs of silicates by the contribution of DBs. The article in itself is ended 
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with a summary. Two appendices give some detail on phonon and breather 
statistics respectively. 

2 From layered silicates to disilicate crystal struc- 
tures 

The synthesis of disilicates from layered silicates, expandable as clay miner- 
als or non-expandable as mica, actually means a reconstructive transforma- 
tion as shown below. 

Layered silicates are made up from two basic building blocks: a sheet of 
edge-sharing [SiOJ units, the tetrahedral sheet, and another one of edge- 
sharing [MOe], the octahedral sheet. There are three main groups of layered 
silicate minerals, according to the combinations of tetrahedral and octahe- 
dral sheets: 1:1, 2:1 and 2:1:1. In 2:1, one octahedral sheet is sandwiched 
between the apices of two tetrahedral ones. In these, also called T-O-T 
silicates, layers are either held together by weak van der Waals forces if 
they are neutral, or may have cations between them for charge balance if 
substitutions in either tetrahedral or octahedral sheet result in a net layer 
charge. The 2:1 layered silicates are classified as trioctahedral or dioctahe- 
dral, after the full occupation of the octahedral sheet by Mg(II) or two thirds 
by Al(III). Talc (trioctahedral) and pyrophillite (dioctahedral) are minerals 
with non-charged layers. In the case of low charge, it results that the clay 
minerals have the capacity to expand by taking up H2O molecules in the 
interlayer space. For high charge, there is mica: phlogopite (trioctahedral) 
and muscovite (dioctahedral). 

Muscovite is a mica which layer charge comes from the isomorphic sub- 
stitution of silicon by aluminium in the tetrahedral sheet. 14 The potassium 
located in the interlayer space, for charge balance, cannot be hydrated; thus, 
muscovite does not expand. Its structure is depicted at the left of Fig. ^ 
At its right, the interlayer space of muscovite is illustrated and it can be 
observed that both surfaces of the upper and lower tetrahedral sheet are 
formed by the basal oxygen atoms from the [SiOJ tetrahedra, which form 
a rough hexagonal honeycomb structure. The interlayer balancing sheet is 
therefore sandwiched in between, potassium occupying the dimples left at 
the centre of each pair of hexagonal cells. In real crystals, although the on- 
site potential created by the silicate layer tends to preserve the symmetry, 
there are always distortions like tetrahedral rotation (the actual situation 
in muscovite is shown in this figure). Each potassium ion is surrounded by 
six other potassium ions in the interlayer sheet with a ditrigonal symmetry. 
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Figure 1: Crystal unit cell of muscovite ICSD 34406. The circles represent 
the potassium ions forming the interlayer sheet, (a = 5.19 A; b = 9.02 A; 
c = 20.0 A; /3=95.7°) 

The distance between the potassium and the basal oxygens layers is 1.45 A. 

Essentially, in both clay minerals and mica [SiOJ tetrahedra are linked 
into infinite two-dimensional networks by sharing three oxygens. How- 
ever, disilicates, or pyrosilicates, are the simplest of the condensed forms, 
where only two tetrahedra share one edge and constitutes the anion Si2 07 6_ 
(Fig|2J). The transformation of any layered-silicate into disilicate involves the 
rupture of two silicon-oxygen bonds by each [SiOJ, whatever the reaction 
mechanism might be, the transformation being reconstructive. 

It is well-known that natural pyrosilicates show a wide range of Si-O-Si 
angle, from 131° to 180°, 15 and that the activation energy is reduced if the 
surface and strain energy terms are diminished by good lattice matching 
across the interface between the new and parent phase. However, this is 
not the present case. It rather seems that the disruption of the tetrahedral 
sheet could be the consequence of localized nonlinear vibration modes as 
commented in the introduction and explained in the following sections. 

3 Experimental RE-disilicates synthesis 

The method used by us to synthesize RE-disilicates consists of a hydrother- 
mal reconstructive process at low temperatures in an isolated reaction vessel 
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Figure 2: Structure of lutetium disilicate. The circles represent the lutetium 
ions. 



constructed in our laboratory. A layered silicate and an aqueous solution are 
the silicon and RE(III) sources, respectively. Up to now, a set of expandable 
clay minerals had been studied, rendering conclusions on the relationship 
between mineralogical compositions and reactivity. 3 The reaction tempera- 
tures were always below the critical one of water; thus both vapour and liq- 
uid phases coexist throughout the whole reaction. Reconstructive structural 
changes occurring in the layered silicate are always analyzed studying the 
long-range order by X-ray powder diffraction (XRD), the chemical environ- 
ment of the main constituent elements of the lattice by magic-angle spinning 
nuclear magnetic resonance spectroscopy (MAS-NMR) and the microstruc- 
tural and microchemical composition by electron microscopy (SEM) and by 
energy-dispersive X-ray (EDX). 

As it has been already mentioned, it has been suggested that DBs might 
bring about an apparent violation of Arrhenius law, leading to chemical re- 
actions being observable at much lower temperatures than expected. 3 It 
has also been suggested that DBs in the interlayer potassium sheet may be 
responsible for the dark lines observed in crystals of mica muscovite 12 ' 16 and 
the existence of moving breathers in mica has been recently proven through 
a sputtering experiment. 13 It implies a LTRT process occurring in a non- 
expendable layered silicate, which contradicts the mechanism published by 
some of the present authors to explain the synthesis of RE-disilicates as as- 
sociated to the capacity for expanding. 3 In order to confirm experimentally 
the hypothesis that the LTRT phenomenon is not an exclusive feature of clay 
minerals, but can also be attributed to non-expandable layered silicates, we 
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have performed the hydrothermal synthesis of Lu-disilicate from muscovite, 
under the same experimental conditions as those used for clay minerals. 

When muscovite is hydrothermally treated in stainless steel reactors at 
300 °C for 72 hours with lutetium nitrate 0.05 M solution, it gives rise to 
Lu2Si20y. Figure^ shows the XRD diagram of the untreated muscovite. It 
reveals numerous hkl basal reflections compatible with the 2Mi polytype and 
a perfect ordering of the layers. After the hydrothermal treatment, the XRD 
pattern shown in figure exhibits a number of specific reflections which 
are consistent with the development of a new crystalline phase Lu2Si2C>7 
(JCPDS file number 76-1871). The SEM photograph of the mica shows big 
flakes whereas the sample submitted to hydrothermal treatment, in addition 
reveals irregular, rough particles, which correspond to Lu2Si2C>7. The 29 Si 
MAS NMR spectra of both expandable and non-expandable layered sili- 
cates show an evolution from Q 3 silicon environment to Q 1 environment. 17 
This result supports that the LTRT phenomenon is common to expandable 
and non-expandable layered silicates, both containing an interlayer sheet of 
cations to balance the layer charge. 

Note that although the hydrothermal treatment of layered silicates, mus- 
covite and others, is, effectively, at the origin of the development of a new 
crystalline phase, as a necessary experimental condition, it is not sufficient. 
In the case of a silicon source different from a layered silicate, which has 
been used by the authors for the first time, the process does not occur. An 
example is illustrated in the Fig. 0J This figure shows the 29 Si MAS NMR 
of SiC>2 submitted to a hydrothermal treatment at 300 °C (b) and the typ- 
ical spectrum of the lutetium disilicate (a). The signal of the SiC>2 after 
hydrothermal treatment remains in the chemical shift range typical of Q 4 
environment (typical of Si02), which is clearly different for Q 1 environment 
(typical for the new phase Lu2Si207). Thus, it demonstrates that under 
hydrothermal condition, SiC>2 does not transform in a new phase. 

4 The conventional chemical kinetics approach 

Transformation processes of minerals in which there is a major reorganiza- 
tion with bonds and, even, change in the chemical composition are classified 
as reconstructive. These obey mechanisms which involve very high activa- 
tion energies when the rupture of strong bonds are involved. Synthesis of 
RE-disilicates from layered silicates requires the rupture of silicon-oxygen 
bonds, which are considerably stronger than the bonds between any other 
element and oxygen. Silicate minerals make up the vast majority of rocks 
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2 THETA 



Figure 3: XRD pattern and SEM micrography of untreated (a) and hy- 
drothermally treated (b) muscovite. m=muscovite, *=Lu2Si207. The com- 
position of the rough particle inside the circle is compatible with Lu2Si2C>7 

and their reconstructive transformation processes show activation energies 
as high as 200 kJ/mol or even higher. 18 Therefore, these transformations 
can be observed in silicate-based minerals at temperatures higher than ap- 
proximately 1000 °C and are apparently impossible at lower ones. 

It is well known that the over-all effect on temperature on the reaction 
rate constant k is expressed by Arrhenius law: 

k = A exp(-EjRT) (1) 

where A and E a are the frequency factor and the activation energy, respec- 
tively. According to it, and without any further insight into the mechanism, 
the rate of a reaction is about 10 9 times faster at 1000 °C than it is at 300 °C 
for an activation energy of 200 kJ/mol. 

It can be shown that the second parameter in the Arrhenius law, namely, 
the frequency factor A, does not support LTRT synthesis of RE-disilicates 
as well. Usually, reaction rates are described within the frame of the quasi- 
equilibrium activated state model. 19 According to the latter, the necessary 
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condition for a transformation to take place at a measurable rate is that 
a sufficient number of atoms have enough energy to achieve the transition 
state. This energy is supplied by thermal fluctuations. The rate of reaction 
is then simply the number of activated complexes passing per second over 
the potential barrier. 

Applying the conventional transition state (activated complex) theory 20-2 
to our case, the simplest formulation of the mechanism can be cast in the 
following form: 

S (layered silicate) +RE(III)(aq) «-> [S-RE(III)]* 

[S-RE(III)]* -> RE 2 Si 2 7 (2) 

The rate constant k for the reaction can be derived by assuming that the 
transition state (or activated complex) is in equilibrium with the reactants. 
If C* represents the concentration of the transition state then the equilibrium 
constant is: 

K* = F » (3) 

[RE(III)] v ' 

The rate constant k and the equilibrium constant K* are related by the 
Eyring equation: 

k = vK*, (4) 

with v = kg T/h, h being the Planck constant. 

To be precise, the rhs of the above expression should be multiplied by a 
factor rj, the transmission coefficient, which is the probability that the com- 
plex will dissociate into products instead of back into reactants. For most 
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reactions rj is between 0.5 and 1.0. Through a thermodynamic formulation 
of K*, it results: 



k = u V ex.p(-AG°*/RT) = 

v V exp(AS */R)exp{-AH°*/RT) = Aexp(—E a /RT) , (5) 

where the superindex stands for normal conditions. In liquid and solid 
systems, the pAV°* term is negligible and AH = AE°* = E a . 

At 300° C we have calculated the factor A by substitution of all the 
parameters in the above expression and it takes the usual value 10 13 — 10 14 
s _1 for a first order reaction. It does not explain the observation of LTRT 
synthesis of RE-disilicates as previously concluded from E a . 

It is well known that, in parallel to the reorganization of the clay, there 
must be nucleation of RE2Si20 7 crystals and that the activation energy 
is reduced if the surface and strain energy terms are diminished by good 
lattice matching across the interface between the new and parent phase. 4 
However, this is not the present case. It rather seems that the disruption 
of the tetrahedral sheet is the consequence of localized nonlinear vibration 
modes, as suggested in Ref. 7. If the vibration modes were delocalized, the 
relationship shown in Ref. 18 between bond angle, Si-0 distance and free 
energy would be incompatible with an appreciable parent structure-directing 
character. 

Ytrium disilicate (Y2Si20 7 ) has four polymorphs, namely y, (3, 7 and 
5. Our LTRT synthesis from the two layered silicates saponite and laponite 
have shown a structure-directing character of the parent clay by only giving 
y-Y2Si207 and <5-Y2Si20 7 , the lower and higher temperature polymorphs. 4 
The relative position of the two tetrahedra in the disilicate unit of the y- 
and 7- polymorphs are similar to their position in the tetrahedral sheet 
of saponite and laponite, which is ~ 141°. The Si-O-Si bond angles in 
Y 2 Si 2 7 polymorphs are 134° (y), 180° (/3), 170° (7) and 158° (8). The 
percentages of y- and 5- phases in the cases of saponite and laponite are 
explained in Ref. 4 as related with the presence of Al(III) in the precursor 
framework. Moreover, the importance of maintaining the local Si-O-Si bond 
angle of precursor structure is shown by the fact that the ^-polymorph has 
been synthesized at more than 365° C below the stability range shown in the 
phase diagram: 



y-Y 2 Si 2 7 /3-Y 2 Si 2 7 
/3-Y 2 Si 2 7 7-Y 2 Si 2 7 
/?-Y 2 Si 2 7 -► 5-Y 2 Si 2 7 



1050(50) °C 
1350(50) °C 
1500(50) °C 
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At the present stage of DBs application to understanding the LTRT 
phenomenon, other experimental results already published by our group, 
such as the influence of isomorphic substitution in clay minerals on the 
reactivity of the latter, are still awaiting their explanation. We hope that 
the further work on DB in solid state physics will permit to interpret all 
those experimental results in the near future. 

5 Discrete breathers model for potassium vibra- 
tions 

As commented above, it has been predicted that discrete breathers could 
bring about an increase of the reaction rate. 7 They are localized vibrational 
modes in networks of nonlinear oscillators, which have been thoroughly stud- 
ied in the last years. 26-28 They have been observed in experiments 8-11 and 
are thought to play an important role in DNA denaturation. 29 They were 
also suggested to be responsible for dark lines in crystals of muscovite 12,16 
and their existence demonstrated in muscovite through a sputtering exper- 
iment. 13 

In this paper we have considered breathers for the out-of-plane move- 
ments of the potassium cations, i.e. the movements transversal to the cation 
layer. This layer is considered to be a 2D triangular lattice, as Fig. [5] shows. 
The Hamiltonian is given by: 

H = ^ ( 2 + V ( U n) + \k ^ ( U ™ ~ U ™') 2 ) ' ( 6 ) 

n \ n'eNN J 

where m = 39.1 amu is the mass of a potassium cation, k is the elastic 
constant of the cation-cation bond, V(itn) is an on-site potential, and the 
second sum is extended to the nearest-neighbours, as indicated in Fig.[5J The 
value of k, i.e. the elastic constant for potassium-potassium bond stretching 
is 10 ± 1 N/m, after Ref. 30. 

A band at 143 cm -1 is identified in Ref. 31 with the K + vibration perpen- 
dicular to the K + -plane in spectra from 30 to 230 cm -1 . We have performed 
infrared spectra in CNRS-LADIR 1 above 200 cm . There appear bands 
at 260, 350 and 420 cm -1 which we assign tentatively to higher order transi- 
tions of the same vibration. Using standard numerical methods to solve the 
Schrodinger equation for the K + vibrations with a potential composed of 
the linear combination with three gaussians and a polynomial of degree six, 

1 Laboratoire de Dynamique, Interactions et Reactivate at CNRS-Thiais, Paris 
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(n,n-1) (n,n) (n,n+1) 



(n+1,n-1) (n+1,n) 



Figure 5: Scheme of the 2D lattice used for discrete breather calculation, ai 
and a2 are the base vectors. 



with adjustable parameters, we have been able to find a suitable potential 
that fits these bands and their intensities. It is given by: 



with D = 453.11 cm" 1 , b 2 = 36.0023 A~ 2 and 7 = 49884 cm- 1 !- 6 . To 
completely define the potential, the limitation to the K + displacement due 
to the muscovite structure has also been taken into account. Fig. shows 
the adjusted potential together with the observed bands. 

The nonlinearity of the on-site potential allows the existence of intrin- 
sic localized modes or discrete breathers apart from low-amplitude linear 
modes. 7 These localized solutions exist as long as no integer multiple of 
their frequency resonates with the frequency of a linear mode and can be 
obtained numerically with machine precision. They are calculated using 
procedures based on the anti-continuous limit (n = 0), 32 which consists of 
calculating an orbit of an isolated oscillator submitted to the potential V(u) 
with a fixed frequency u^, and using this solution as a seed for calculating 
the solution of the full dynamical equations (k 7^ 0) through a continuation 
method. 

In principle, we have considered excitations in the soft part of the on-site 
potential (the energy of an isolated oscillator decreases with the frequency) 




(7) 
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-1 -0.5 0.5 1 

x (Angstrom) 



Figure 6: Numerically calculated on-site potential for the potassium vibra- 
tion based in the observed bands (horizontal lines) 



because this is the region that fits the experimental data. In this case, the 
frequency of a breather is below the linear modes band and none of its 
harmonics coincide with the frequency of a linear mode for a given value of 

K. 

The linear modes are solutions of the equations: 

mu n}n > +mU%U n , n t - K (Un t n'-1 + U n ,n'+l + «n-l,n'-l + «n-l,n' + 

Un+l,n'-l + U n +l,n' ~ 6Un,n') = , (8) 

with lo 2 = V"(0)/m. Then, the linear modes spectrum is the following: 

oo 2 = lo 2 [1 + AC (sin 2 (^/2) + sin 2 (0 2 /2) + sin 2 (0 2 /2 - ^/2))], (9) 

with C = K/(mu> 2 ) 0.15 and d{ G [— vr,7r]. Thus, the frequency of the 
linear modes lies in the interval v € ^n(l, Vl + 9 C), with v§ = u /(2tt) = 
167.50 cm -1 . In consequence, no multiple of must lie in this band. 

Fig. Q shows the profiles of two breathers with different frequencies and 
Fig. |H] shows the dependence of the breather energy with respect to the 
frequency. We can observe that there is a minimum energy A ~ 23 kJ/ mol 
for breather creation as it is known to happen for 2-D and 3-D breathers. 33 
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These breathers are single symmetric breathers, but there are many 
other different types, as breathers with different symmetries, multibreathers, 
breathers above the phonon band, which will have different energies, range 
of existence and probability of excitation. For example, symmetric single 
breathers with frequencies above the phonon band have energies between 
240 and 500 kJ/mol. The obtention of a complete description for all the 
different breather types and their ranges of existence and stability is a long 
task that we do not pursue here. 




Figure 7: (Left) Energy density profile for a breather with frequency v\> = 
0.97^o an d energy 25.6 kJ/mol. (Right) The same but for a breather with 
frequency fb = 0.85 uq and energy 36.3 kJ/mol. Notice that the localization 
is higher for the lower frequency. The x and y coordinates are in lattice 
units of 5.2 A, vq = 167.5 cm -1 . 



6 Breather statistics and its effect on the reaction 
rate 

6.1 Breather statistics theory 

The effect of the temperature on the reaction rate constant k is given by the 
Arrhenius law k = A exp(— /3 E a ), with f3 = 1/RT. The Boltzmann factor 
exp(— (3 E a ) represents the fraction of vibrational modes that are able to 
deliver the activation energy E a for the reaction to proceed. The frequency 
factor A depends on the reaction kinetics and although it is difficult to 
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Figure 8: Dependence of the breather energies with the relative frequency 
for symmetric single breathers with frequencies below the phonon band. The 
dashed line indicates that breathers are unstable, uq = 167.5 cm -1 . 



calculate the estimation given above is 10 — 10 s , or the order of 
magnitude of most chemical reactions. 

The statistics of breathers is not at all completed and a full statistical 
analysis of muscovite and the reaction kinetics is a formidable task. Here 
we use the theory developed for 2-D breathers by Piazza et afi 4 which is 
described in the first subsection of Appendix El with some modifications in 
subsection IB.2I Note that this theory is not deduced from first principles 
but it is based in heuristic reasoning. It has to be considered a reasonable 
approximation justified by its good accuracy with numerical simulations. 
Here we summarize the main results that are useful to us: 

a) The probability density P\ } (E), defined so as Pb{E) dE is the proba- 
bility of existence or the mean fraction of breathers with energies between 
E and E + dE, is given by 

P h (E) dE = J^fiiE - AY exp[-(3(E - A)] , (10) 

with T(z + 1), the gamma function, and z, a parameter which does not 
depend on the energy. The cumulative probability, C^(E) or mean fraction 
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of breathers with energy above or equal to E is given by: 

C h (E) = T(z + l)" 1 7 (z + 1,(3(E- A)) , (11) 

where j(z, e) is the incomplete gamma function, defined as y 2_1 exp(— y) dy. 

b) A population of large breathers tends to develop with mean energy 
(E) = A + (z + 1) IcbT and maximum probability for E(P raax ) = A + z k^T. 

Fig. shows the probabilities and cumulative probabilities of phonons 
and breathers, the latter for A = 20 kJ/mol and z = 2. Note that for large 
energies of the order of the activation energies, C^(E), although small, is 
several orders of magnitude larger than C p h(E). The drawback is that the 
number of breathers per site is much smaller than the number of phonons, 
with typical values of 1CP 3 obtained in numerical simulations. 

This theory can be modified to take into account that some breather 
types have an upper limit for their energies, due to bifurcations or instabili- 
ties as explained in IB.ll The modified probability densities and cumulative 
probabilities are given in Eqs. (|B-6IB-7|) . 

6.2 Numerical simulations 

We have performed a number of numerical simulations with the following 
procedure: a) A network of 50 x 50 oscillators is given a number of times 
(about 500) the same energy with different random distributions of veloci- 
ties; b) They are left evolve a sufficiently long time so as they attain states of 
thermal equilibrium, for which the mean temperature for all the simulations 
is calculated; c) The networks are cooled by adding a dissipative term to the 
dynamical equations for the oscillators at its borders until the energy stops 
diminishing; d) The number of breathers is counted and their energies are 
calculated. The number of breathers divided by the number of sites and the 
number of simulations gives the mean number of breathers (n^). 

Note that the temperature can not be fixed from the outcome. For some 
energy we obtain a temperature Tc = 280.25 C, approximately the experi- 
mental one, and the data given here correspond to that case. We obtain a 
mean number of breathers per site (n h ) = 0.92 x 10~ 3 , with a mean energy 
(E) ~ 70 kJ/mol which would correspond to z = 9.19 at the experimental 
temperature. However, the probability density Pb(E) does not correspond 
to a curve as the one given by Eq. ([TO]) and plotted in Figure ITU because 
it is a curve with several maxima and minima (see Fig. |SJ). The reason is 
that with the numerical simulations, we obtain all types of multibreathers 
and single breathers with different symmetries and not only the symmetric, 
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single, exact breathers obtained numerically in the previous section, which 
are continuation from a single excited oscillator at the anticontinuous limit. 

Therefore, the numerical curve Pb(E) can be seen as a numerical spec- 
trum for the different breather energies and forms of vibration. The ideal 
objective, as with other types of spectra, would be to know each type of 
breather, with its dispersion curve E(i%), and the relative probability of its 
appearance and to be able to reproduce exactly the numerical spectrum. 

In principle, all the breather types could be obtained exactly with dif- 
ferent conditions at the anticontinuous limit with different frequencies and 
by path continuation by changing the frequency and studying the different 
branches at the possible bifurcations. This is a long and difficult task that 
we do not pursue here. Instead, we try to fit approximately the numerical 
Pb(E), with a small number of breather types, each one characterized by its 
minimum energy A, maximum energy Eyi, which can be do, and parameter 
z (see IB. H and IB.2| . each breather type with a different probability to occur. 
In this way, we know that we cannot fit exactly the numerical spectrum be- 
cause we are most probably substituting a number of breather types with an 
average one. In any case, our numerical Pb(E) is also an approximation, as 
we would need a extremely large number of simulations to obtain the actual 
curve. 

Note that the breather spectrum will not appear in an experimental one 
for three reasons: a) The number of breathers is about 10 -3 the number 
of phonons, and the spectrum is basically dominated by the one-phonon 
transitions; b) Breathers are localized and, therefore, they cannot be excited 
by infrared, raman of neutron spectra. 36 

Figures |H1 and ^] show the numerical and theoretical, probability den- 
sities and cumulative probabilities, respectively. The parameters of the 
breathers are: 

A (kJ/mol) 23.9 36.6 41.4 62.2 67.3 82.9 
z 1.50 1.17 3.00 0.52 2.07 1.80 

E M (kJ/mol)) 46.9 94.4 

Probability 0.103 0.026 0.281 0.097 0.202 0.290 

6.3 Effect on the reaction rate 

According to Ref. 18 the lowest estimates of the activation energy for a 
reconstructive transformation as the one described above are E & = 100 — 
200kJ/mol. Let n p h(£'a) — exp(— 0E a ) and n^Ea) be the mean number 
of phonons and breathers, respectively, per site, with energies E > E a . 
Then, n h (E a ) = {n h )C h (E a ), with (n h ) ~ 0.92 • 10" 3 , the mean number of 
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Figure 9: Breathers spectra, i.e., breather probability densities Pb(E) ob- 
tained numerically (line with dots) and theoretically (continuous line). The 
latter is obtained by considering six different types of breathers. See text. 



breathers per site obtained numerically, and Cb(E a ), the cumulative proba- 
bility described in the previous subsection, obtained with six different types 
of breathers in order to fit the numerical probability density. The ratio of 
the number of breathers to phonons becomes nb(-E a )/n p h(-E 10 4 - 10 5 . 

The reaction rate constant with breathers would be k\> = A^n^Eg), 
with A\>, the frequency factor for breathers A^, which should be different 
from A, the frequency factor for phonons. We will assume that A^ = A for 
the purpose of comparison. The ratio of reaction rates for breathers and 
phonons would be k h /k = n h (E a )/n ph (E a ) ~ 10 4 - 10 5 for E a = 100 - 200 
kJ/mol. In other words, as the three days experimental time leads to about 
30% of the transformation performed, the time without breathers to obtain 
the same result, would be 10 4 — 10 5 times larger and, thus, completely 
unobservable. 

Two factors are likely to increase further the reaction rate with breathers. 
First, since a discrete breather is strongly localized, it seems much more 
capable of delivering the energy for breaking a Si-0 bond, which implies that 
Ah should in effect be much larger than A. Second, for larger systems than 
the one used in our simulations, the fluctuations may have larger energies 
and, therefore, excite other types of breathers with higher energies, as, for 



17 



1.4 



1.2 




20 40 60 80 100 120 

E (kJ/mol) 



Figure 10: Breather cumulative probabilities C\^{E) obtained numerically 
(line with dots) and theoretically (continuous line). The latter considering 
six different types of breathers. See text. 



example, breathers with frequencies above the phonon band, which have 
energies between 240 and 500 kJ/mol. These breathers would increase the 
fraction of breathers above the activation energy and therefore the reaction 
rate. 

7 Summary and conclusions 

Low temperature reconstructive transformations (LTRT) have been achieved 
in layered silicated by some of the authors at temperatures about 600° C 
lower than previously reported. This is a phenomenon for which there is 
presently no plausible explanation since the bonds involved are the same 
as in other transformations. New experiments performed by some of the 
authors on mica muscovite, a non-expandable silicate have discarded their 
previous hypothesis of LTRT been caused by the expansion of the intersheet 
layer. 

We have constructed a model for breathers in the cation later, for which 
we have obtained reasonable parameters, and with a mixture of numerics and 
theory we have estimated their effect in the reaction rate. The results are 
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that they would increase enormously the reaction rate and, thus, explain the 
observed LTRT. This can be easily explained by the fact that, although there 
are much less breathers than phonons, there are many more with energies 
above the expected activation energy. Certainly, the statistical theory of 
breathers is only an approximation, and the numerics cannot be precise at 
larger energies for which there are so few breathers, except if an enormous 
number of simulations could be performed. However, an established fact in 
breather theory is that large breathers have longer life time than small ones 
and thus tend to overpopulate the regions of high energies if compared to 
Maxwell-Boltzmann statistics for phonons. Moreover, as they are localized, 
it seems that they can deliver more easily the required energy to break bonds. 
The sum of this facts, i.e., localization, much higher number of breathers 
above a given activation energy and apparent diminution of the activation 
energy suggests that DBs are good candidates to explain LTRTs. 
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Appendices 

A Phonon statistics 

A.l Quantum statistics of one oscillator 

A quantum harmonic oscillator with frequency u, equal to its classical one, 
has energies E n = (n + 5) hu. If in contact with a thermal bath at tem- 
perature T, the probability that it has energy E n is given by P(E n ) = 
A exp(- [3 E n ) = Aexp(-/3 (n + \)Ku), with (3 = l/k B T, k B being the 
Boltzmann constant. 

A can be obtained by the normalization condition YlnLoPn = A Z = 1, 
where Z = ex P[ — P (n + ^) hu>] is known as the partition function for 

the oscillator. Z is a geometric series which can be easily summed leading 
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to: 

= FV l-g-i- . (A-1) 

Therefore A = l/Z and 

P n = P(£7 n ) = ex P[~^(^+ 5)M = eX p(-/?reM[l _ exp(-^w)] (A-2) 

Note that n is the excitation number of the oscillator, however, in solid 
state physics, where ui is the frequency of a normal mode of the solid, it 
is customary to speak of n as the number of phonons with frequency to. 
Once known P n a number of quantities can be readily calculated. The 
mean energy (E)=Z- 1 £^° =0 (n + \)hw exp[-/3(n + \)hw] = -Z^dZ/dfi = 
—d\og(Z)/d(3, which leads to: 

(E) = (- + — L- )hu. (A-3) 

X ' V 2 exp(/3^j)-l ; v ' 

Consequently, the mean excitation number or mean number of phonons 

is: 

{n) = ex P(/3 L)-r (A " 4) 

For high temperatures huj/ksT « 1, the mean energy becomes (E) ~ 
keT which is the classical mean energy of a harmonic oscillator and (re) ~ 

Of particular importance for the present problem is the cumulative prob- 
ability C(E a ), i.e., the probability that the oscillator has energy E a or 
higher above the ground state hcu/2, i.e., that it can deliver the energy 
E a . Let m be the minimum integer so as m?iw > E a , i.e., m = \E a /hw~\, 
with \x], the ceiling function that rounds x towards plus infinity. Then 
C(E a ) = Z- 1 E~ m exp[-/3(i+n)M = n° =0 exp[-/3(m + I + n)M 
= Z~ x exp(— fimhw) x J2^=o ex P[ — 0{\ + n)ftw\, therefore: 

C(E a ) = exp(-/?m/^) = exp(-/3[^] hw) . (A-5) 

If m is of the order of a few tens, the expression above approaches to the clas- 
sical expression C c i ass (-E a ) = exp(— (3E a ). Note, however, that the quantum 
probability is somewhat smaller. The ratio between the quantum probabil- 
ity and the classical one is between exp(— flhiS) and 1. For a frequency as 
the one given here for the phonons in the K + plane ujq = 3.16 • 10 13 s~ 1 and 
T = 573K this ratio is exp(-/9/iwo) ~ 0.65 and 0.81 for T = 1173K. 
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A. 2 Normal modes and phonons 

Let us consider a solid in the linear approximation, with Nf degrees of free- 
dom, with Nf = 3 x N a for a three-dimensional solid with N a atoms. There 
are Nf normal modes with frequencies tOi and wave numbers k;, each one 
equivalent to a independent harmonic oscillator, and therefore, the previous 
subsection can be applied to it. The properties of the solid are simply the 
sum of the properties of the isolated oscillators, just adding the subindex i 
to the formulae in the preceding section and summing up, i.e., the excitation 
number (rii) of the mode i (or the number of phonons) and the energy of 
the solid are given by: 

<"•> = ex P (fflL.) + 1 ; B "BW + -^ (A-6) 

1 = 1 

The number of modes N p h(E a ) with energy larger or equal to E a above 
their ground state is 

N ph (E>E a ) = J2exp(-f3\^]hj l ). (A-7) 
7=1 n ^ 

If E a is about a few tens larger than any fUdi, then we simply have: 

N ph (E > E a ) ~ N f exp(-/3£ a ) , (A-8) 

which is the classical expression. Again, the ratio between the quantum and 
classical expressions of N p h(E > E a ) is somewhat smaller than the unity. 

Therefore, the cumulative probability C(E Su ) i.e., the fraction of modes 
with energies greater of equal to E a , becomes: 

C{E a ) ~ exp(-/3£ a ) . (A-9) 

B Breather statistics 

B. l Breathers with hard on-site potential 

The breather statistics theory developed in Ref. 34 for 2D breather in a sys- 
tem with hard on-site potential is based is some simple hypotheses, which, 
in principle, can be fairly general: 

1 . An established fact is that breathers in two and three dimensions have 
a minimum energy A. 33 
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2. The rate of creation of breathers with energy E, B(E), is propor- 
tional to exp(— (3E), since breathers form from fluctuations through 
an activation process. 

3. The probability per unit time that a breather with energy E is de- 
stroyed, D(E), is postulated to be inversely proportional to (E — A) z , 
with z a constant (which means as other constants hereafter that it 
does not change with the energy E) that depends on the system. This 
law is the simplest mathematical expression that takes into account 
that large breathers have longer lives than smaller ones, with the life- 
time of breathers with minimum energy A equal to zero. It has to be 
considered as an approximation as it is not derived from first princi- 
ples. 

Let Pb (E) dE be the probability of existence (or the mean fraction) of 
breathers with energy between E and E + dE. The rate of destruction 
of breathers with energy E is proportional to D(E) and P^(E), therefore, 
exp(-0E) = AP h (E)(E - A)~ z , A being a constant, or, AP h (E)=(E - 
A) z exp(— (3E). Since P b (E)dE = 1, A can be obtained using the change 
of variable y = fi(E - A): 

poo 

A = {E-A f exp(-0E) = 
J A 

6XP /? ( ~ / ? A) jH V z exp(-y)dy = exp(-/3A)/r (z+1) r(z + 1) , (B-l) 

with Y(z + 1) = / °° y z exp(— y)dy, the Gamma function. Thus, the proba- 
bility of breathers with energy between E and E + dE is given by: 



A y ' * K H ' exp(-/3A)/3-(^+ 1 )r(z + l) 



The mean energy is given by: 



■(E — A) z exp[—/3(E — A)] . (B-2) 



f'OG f'OG 

(E) = / EP h (E)dE = A + / (E- A)P b (E)dE 
J A J A 

poo Rz+1 

A + J A YJ^iE-^eM-PiE-AME 
A + — 1 y z + l exp(-y)dE = A + T{z + 2) 



PT(z + l)J * ^ a > /3r(z+l) 

A + (z + l)k B T. (B-3) 
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The cumulative probability C\,(E), i.e., the probability that a breather 
has energy higher than E, is given by: 



/oo poo R z <i 

P b (E)<US = j ^ + l) (E~ A) 2 exp[-/?(£ - A)] dE = 

1 y z exp(-y)dy = -i -r - , (B-4) 



r(* + 1) Jke_ A ) * y ' r(* + i) 

where T(z + l,x) = y z exp(— y)dy is the first incomplete Gamma func- 
tion. 35 

The energy for which the probability is maximum is given by E(P max ) = 
A+z ke T, which shows that breathers tend to populate higher energies than 
phonons. As an example, Fig.^Jshows Pb(E) and Cb(E) for breathers with 
A = 20 kJ/mol and z = 2 and the equivalent magnitudes for phonons. The 
large energies the larger values of Cb(E) soon compensate for the much 
smaller number of breathers than phonons (around 10 -3 ). Although the 
extrapolation of C\>(E) to large energies has to be done with caution and a 
more elaborate theory has yet to be developed, the basic fact that breathers 
tend to populate higher energies than phonons can be accepted. Piazza et 
al have succeeded in fitting the cumulative probability in Eq. lB-4l with the 
observed one in numerical experiments, which proves that the hypotheses 
above are reasonable. 



B.2 Breathers with maximum energy 

Hereafter we modify slightly their theory developed above. The on-site 
potential for the 2D system in Ref. 34 is hard, i.e., the frequency of the 
isolated oscillators increases with the frequency, with the consequences that 
the breather frequency lies above the phonon band, its energy increases with 
its frequency and can be considered as unbounded. 

For soft on-site potentials or, as in the present paper, potentials with 
both soft and hard parts, breather energies may have an upper limit be- 
cause the breather frequency enters the phonon band or because a bifur- 
cation where the breather disappears or transforms into a different one, as 
a multibreather or a breather with different symmetries. The changes are 
obtained by introducing an upper limit Em in the integrals with respect to 
the energy. Then, Eq. (|B-1|) becomes (with y = j3 (E — A)): 

A = (E- A) z exp(-/3E) = > J y* exp(-y)dy = 

exp(/?A)/r (2+1) 7(^ + l,/3 (E M - A)) (B-5) 
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Figure 11: Comparison of the phonon and breather probabilities densities 
(left) and cumulative probabilities (right). The breather values have been 
obtained for z = 2 and A = 20 kJ/mol. The temperature is T = 600 K and 
IcbT ~ 5 kJ/mol. 



where j(z + 1, x) = y z exp(— y) dy is the second incomplete gamma func- 
tion. 35 

Therefore, the probability density becomes: 

P h (E) - -{E - A) eM-PE) l{z + h(3[EM _ A)) • ( B -6) 

The cumulative probability becomes: 

M ' Je K 1 Je j(z+1,(3(E m -A)) 

For Em >> A, the expressions above for Pb(E) and Cb(E) transform into 
the expressions in Eqs. (|B-.3|B-4jl 
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In a system like ours there are different types of breathers and the prob- 
abilities or cumulative probabilities calculated above correspond to each 
type with different minimum and maximum energies A and Em (or without 
maximum energy), and parameter z. The total number of breathers and the 
relative probability of each type are unsolved questions. The latter probably 
depends on the temperature, the breather energies, the phase space occu- 
pied by each breather and its equivalent ones through symmetries, and the 
breather profile, which might be excited more or less easily by the phonons. 
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